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1 Introduction 

Electronic transport in semiconductor nanostructures is one of the central issues of research in modern 
condensed matter physics, see, e.g., [UEl for reviews. In this article, we review recent results on transport 
of two-dimensional electron gases (2DEG) in moderately strong transverse magnetic fields B. Specifically, 
we concentrate on a range of B which are classically strong (i.e., w c T tr ^> 1, where lo c is the cyclotron 
frequency and r tr the transport relaxation time), but where quantum localization effects (and, correspond- 
ingly, quantum Hall physics) are not developed yet. There exists a broad class of phenomena that lead to 
strong magnetoresistivity p xx (B) in this range of fields, in view of the developed cyclotron motion. These 
phenomena are discussed in the present review. 

2 Quasiclassical memory effects in magnetoresistance 

The recent interest in quasiclassical transport properties of a 2DEG has been largely motivated by the 
experimental and practical importance of high-mobility heterostructures, in which charged impurities are 
separated from the 2DEG by a wide spacer. The correlation radius of disorder produced by the impurities is 
usually much larger than the Fermi wave length of electrons and transport in the 2DEG retains signatures 
of the underlying quasiclassical dynamics of the particles. On the theoretical side, much of the interest 
has been inspired by a variety of anomalous transport phenomena which, while being essentially classical, 
cannot be described by Boltzmann-Drude kinetic theory. The quasiclassical "non-Boltzmann" phenomena 
in disordered electron systems are due to correlations of scattering acts at the points where quasiclassical 
paths self-intersect, which gives rise to memory effects, neglected in the conventional Boltzmann equation. 
In particular, the non-Markovian kinetics yields a strong magnetoresistance (MR) and anomalies in the ac 
response. The strength of the quasiclassical anomalies depends on the ratio d/l, where d is the correlation 
radius of disorder, I the mean free path, and grows with increasing d as a power of this parameter. Since 
quantum corrections (weak localization, Altshuler-Aronov corrections, etc.) are governed by a different 
small parameter 1 / UfI 1, where kp is the Fermi wave vector, it is the long-range correlations of disorder 
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with kpd ^> 1 that reveal the quasiclassical anomalies. In this section, we focus on the quasiclassical 
memory effects in magnetotransport. 

2. 1 Magnetoresistance of a 2DEG subject to smooth disorder 

We begin by considering the MR of a 2DEG in the presence of Gaussian disorder which is smooth on the 
scale of kp 1 . First of all, we recall that there exists a finite MR OH) even within the collision-integral 
approximation, the source of which is the bending of quasiclassical trajectories by magnetic field B on the 
scale of the correlation radius d. This leads to a small negative MR, Ap xx /po ~ — (d/R c ) 2 , where po is 
the Drude resistivity and R c is the cyclotron radius. Remarkably, the non-Markovian kinetics gives rise to 
a much stronger positive MR 0, which may even be much larger than unity. 

To systematically treat the quasiclassical memory effects, the starting point is the disorder-averaged 
expression for the conductivity tensor a in terms of the exact Liouville operator L: 




Here v is the density of states, vp the Fermi velocity, 4> the velocity angle on the Fermi surface. The 
operator L = Lq + SL is given by the sum of the free part Lq = —ioj + v F nW + uj c d^, where uj c is the 
cyclotron frequency, n = (cos </>, sin </>), and the part induced by a random scalar potential V(r), 

SL = Sv{r)nV + [VSv(r)](z x n)<9 , (2) 

where Sv(r) = v(r)—v F denotes a fluctuation of the local Fermi velocity v(r) = [v F — 2V(r)/m)] 1 ^ 2 and 
z is a unit vector in z direction. Expanding Eq. ([TJ in SL, averaging over the disorder and resumming the 
series, the diagonal resistivity p xx — 2(—iuj + M xx )/e 2 vv 2 F is represented in terms of the disorder-induced 
self-energy ("memory function") M xx . 

The Drude result Mix — t^ 1 , where r tr is the momentum relaxation time, follows as the first term in a 
perturbative expansion of the self-energy in the strength of disorder, M X °J = — (SL L^ 1 SL). Substituting 
the propagator L^ 1 renormalized by impurity scattering for L^ 1 in the latter expression yields the main 
contribution to the MR associated with the quasiclassical memory effects. When calculating L^ 1 in the 
case of long-range disorder, the stochastic motion of particles can be approximated by a Fokker-Planck 
equation corresponding to the diffusion in momentum space, so that Ld is written as 

L D = -iw + v F nS7 + u c d^ - r^d 2 ^ . (3) 

The leading correction to M xx due to the self-intersection of quasiclassical paths then reads: 

AM XX = (4ir 3 m 2 v 2 F y 1 fd 2 qdcj) sincf) sin(> - (j) q )q 2 W(q) g D (uj, q, 0) , (4) 

where gn{u}, q, <j>) is the Fourier-transformed real-space solution of the equation £d.9d = sin0sin(</> — 
4> q )S(r), 4> q is the angle of q, and W(q) is the Fourier transform of the correlator (V(0)V(r)). For the 
case of impurities separated from the 2DEG by a spacer of width d, Eq. gives 

Ap xx /p a = AM XX /M^ = 2^-^(3/2) [d/lf (tu c T n -) 9 / 2 . (5) 

One sees that the MR due to the memory effects is much larger than that related to the effect of magnetic 
field on the collision integral for w c r tl ^> (l/d) 2 ^ 5 . The MR (0 becomes of order unity when the mean- 
square shift of the guiding center of a cyclotron orbit after one revolution, 

5 = 2Tr 1 / 2 v F T tl -/(Lo c T tI ) 3 / 2 , (6) 
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becomes of order d, which happens at w c r tl ~ (I /d) 2 ^ 3 . At higher fields, the strong positive MR is followed 
by a sharp (exponential) falloff of p xx with growing B |6|: 



(7) 



which is due to the increasing adiabaticity of the electron dynamics and the related quasiclassical localiza- 
tion. The self-intersection induced MR, given by Eq. ©, may be considered as a precursor of the adiabatic 
localization. 

In the limit of weak inhomogeneities, the return-induced MR depends in an essential way on the behav- 
ior of the disorder under time reversal [5 ]; in particular, it is strongly enhanced in a random magnetic field 
(RMF). The case of a smoothly varying RMF is of particular interest in view of the composite-fermion 
description of the transport properties of a half-filled Landau level [7]. Also, a long-range RMF has been 
realized in semiconductor heterostructures by attaching superconducting or ferromagnetic overlayers or by 
"prepatterning" the sample (randomly curving the 2DEG layer). Following the same route as for the case 
of a random scalar potential, the MR due to the quasiclassical memory effects is obtained as [5] 



pxm/po = 1/2 + [1/4 + {B/B ) 2 ] 1/2 



(8) 



where B is the characteristic amplitude of fluctuations of the RMF. For the composite-fermion model at 
half-filling, (B/B ) 2 is represented as 2(d/l)(ui c T a ) 2 . The adiabatic localization in the RMF begins at 
B ~ B (l/d)^ 6 00, so that there is a wide range of B in which the positive MR ([8]) is strong. 





B/B„ 



Fig. 1 Magnetoresistivity in a random potential from nu- 
merical simulations in comparison with Eq. (J5} for l/d = 
290. 



Fig. 2 Magnetoresistivity in a random magnetic field 
from numerical simulations for three different strengths 
of the disorder a — (eBo / 'mc)(d / 'v f); the full line cor- 
responds to Eq. l[8j. 



The numerically calculated MR [5| for both types of disorder shown in Figs. Q] and [2] confirm the 
theoretically predicted positive MR. Note that the MR in the RMF at moderately small d/l still exists, but 
becomes weak; this is the region of d/l relevant to the composite-fermion model. The numerical data for 
d/l ~ 0.1 — 0.2 agree well [8] with the experimental results [9| for the MR around half-filling. Recent 
unpublished numerical simulations ITOll also confirm the u 9 J 2 positive MR in a smooth potential; however, 
the numerical coefficient is found [ 1 ] to be smaller by a factor ~ 2 as compared to Eq. ([5]). 



2.2 Magnetoresistance of a 2DEG subject to two-component disorder 

We now turn to a 2DEG moving in a non-Gaussian random potential represented by rare strong short- 
range scatterers and subject additionally to a smooth random potential discussed in Sec. 12.11 One of the 
most relevant experimental realizations of the model of "two-component disorder" is random antidot (AD) 
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arrays, where the potential barriers around the ADs can be modeled as hard disks reflecting electrons 
specularly (for experimental work on the dc MR in random AD arrays see, e.g., Refs. ifTTl [121 [131 [141 [151 
[l6lD . The model is also applicable to the description of the MR in an unstructured ultra-high mobility 2DEG 
with a wide spacer, where large-angle scattering on residual interface impurities and interface roughness 
becomes important lfT7l[T8l[T9l . limiting the mobility with further increasing width of the spacer. From the 
theoretical point of view, the interplay of the two types of inhomogeneities is quite remarkable in that it 
yields nontrivial physics which is absent in the limiting cases, when only one type is present. In particular, 
although in the extreme of strong B — > oo the resistivity p xx tends to zero in either of the limiting cases, 
it diverges in the presence of both types of disorder [20|. Also, in the experimentally relevant situation 
when the mean free path at zero B is determined by scattering on ADs, the presence of weak long-range 
disorder will nonetheless become of crucial importance with increasing B [21 20]. The magnetotransport 
in the Lorentz-gas model describing an AD array without smooth disorder has been studied, in particular, 
in Refs. Il22l l23l l24l l25l l26l . A strong negative MR followed by a metal-insulator phase transition in a 
strong magnetic field was found in [22 1 23 25 1; the low-field anomalous MR due to the "corridor effect" 
was investigated in [26|. Smooth disorder, however, changes the MR qualitatively, as discussed below. 

Generalizing the formalism described in Sec. 12. II to the case of two-scale disorder, we represent the 
Liouville operator as L = Ld + 5L, where Ld is given by Eq. (0 an d includes interaction with the smooth 
disorder (r sm is the corresponding transport scattering time), whereas SL = — J^. /r^ describes collisions 
with ADs whose random positions are . The Fourier transform I q of the collision operator /r^ yields the 
transport time rg for the scattering by the AD array of density ns through inn = — n / n ST~s- We assume 
that lo~ 1 < tj < T sm , so that the total transport rate is determined by ADs, t,^ 1 = r^ 1 + r^ 1 ~ tJ . 
The leading contribution to AM XX reads 

AM XX = -n s J (d(j>/ 7r) cos <p IrDI r cos </> , (9) 

where the propagator D includes the first-order self-energy: D — (Ld — n-slo) ■ As compared to the 
Lorentz model 11271 in which only hard-disk scatterers are present, new physics emerges in the limit S ^> a, 
where a is the radius of the ADs and S is the shift of the cyclotron orbit after one revolution due to scattering 
on smooth disorder (cf. Sec. 12. 11 1. In particular, Eq. © yields ETI 

A/WPo = -k/^o) 2 , Lo = {27m s ) 1 / 2 v F (3T S /2T sm ) 1 / 4 , (10) 

for io c <C coo. The mechanism of the negative MR ( fTOb can be understood as follows. If one associates with 
a particle trajectory a strip of width 2a, the ratio (ui c /uio) 2 gives the fraction of the area "explored" twice, 
which implies an effective reduction of the exploration rate and thus a longer time between collisions with 
different ADs. The negative MR ( [Tol l should be contrasted with the positive MR (0 for one-scale smooth 
disorder, where the passages through the same area lead to an enhanced scattering rate. 

For oj c ^> luq, the renormalized scattering time t' s 3> t$ should be found self-consistently from the 
condition ns£,R c ~ 1, where £ ~ 5{vft' s / Rc) 1 ^ 2 is a characteristic end-to-end size of the diffusive 
guiding-center trajectory in time t' s , which gives the B~ 4 falloff [21] with increasing B: 

Pxx/po ~ t s /t' s ~ (T S /T sm )(n s R 2 c ) 2 ■ (11) 

Equation (fTTT i is valid as long as t' s <§; r sm , which is rewritten as ngR 2 ^S> 1. In the opposite limit (but still 
for 8 <C d), the scattering on ADs stops playing any role and p xx has a plateau with 

Pxx/po = Ts/r sm . (12) 

Let us now briefly outline what happens at larger B, namely for i/d < 1 1201 . In dilute AD arrays 
for intermediate B, the exponential slowdown of the electron dynamics induced by adiabatic localization 
transforms into p xx /po ~ nsR c d\n(l/nsR c d) oc B^ 1 \nB due to rare collisions with ADs which mix 
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otherwise closed drift trajectories in phase space. Most interesting, however, is the behavior of p xx in 
the limit of large B, where p xx starts to grow as a power law with increasing B. This behavior can be 
most clearly seen in a "hydrodynamic model" of the chaotic AD array (ns — > oo, t$ — const), where 
the problem can be mapped onto that of advection-diffusion transport 11281 . In the limit B — ► oo the 
hydrodynamic model predicts p xx /po ~ (TgVFd/T sm R^.) 5 ^ 13 oc B 10 / 13 . The physics of the positive MR 
is a percolation of drifting cyclotron orbits limited by scattering on ADs. The growth is checked on the side 
of large B by quantum effects (Shubnikov-de Haas oscillations). The different types of the MR in the two- 
component model 1 20 1 are illustrated in Fig.|3]for various concentrations of short-range scatterers (ADs): 
(i) the MR is positive owing to the "diffusion-controlled percolation"; (ii) due to the adiabatic localization, 
the concentration of conducting electrons decreases as B^ 1 In B before the percolation becomes effective, 
which yields a negative MR p xx (B) oc In B for intermediate B; (iii) an exponentially sharp falloff of 
Pxx(B) at B ~ B a d (shown as a vertical jump) separates the diffusive and drift regimes; (iv) because of 
the memory effects, the collision time for scattering by ADs is increased as compared to the Drude value 
already in the diffusive regime (B -C B a d), which leads to the negative MR p xx (B) oc B~ 4 for small B; 
(v) for intermediate B, the scattering on ADs stops playing any role and p xx (B) is saturated at a value 
determined by the long-range disorder only, whereas at larger fields the diffusion-controlled percolation 
gives rise to a positive MR. 



' Pxx< B > 


(I) 

s£ (H) 




j/^/ (Hi) 




\ {IV) 







Bad B 



Fig. 3 Schematic behavior of the magnetoresistivity 
pxx (B) on a log-log scale in the two-component disorder 
model for different values of the concentration of anti- 
dots n: n w > n (ll) > . . . > n {v \ keeping all other pa- 
rameters (ts, i~sm, d) fixed. Only one characteristic field 
Bad is shown, at which the crossover between diffusive 
dynamics and adiabatic drift in the long-range potential 
takes place. 
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Fig. 4 Magnetoresistivity in a two-component-disorder 
model at fixed ts and different r sm from numerical simu- 
lations; T sm /Ts = oo (Lorentz gas, A), 111 (OX 70 (□), 
37 (0). Inset: uio determined from the fit to the quadratic 
behavior given by Eq. dlOt ; the full line corresponds to 
the analytical result for ujq from Eq. dlOt . 



The MR obtained numerically ETl for the model of two-scale disorder in the regime of Eqs. ( TTOb - 
<n~2l >. Fig. |4] shows good agreement with the above analytical results. As far as the experimental data are 
concerned, for typical parameters [ 13] of AD arrays, ns = (0.6 /im)' 2 , vpTs = 1.3 pm, VFT sm = 16 pm, 
and the electron density 5 x 10 11 cm~ 2 , the field Bo corresponding to the frequency uio given by Eq. ( [Tol l 
is ~ 0.3 T, in agreement with the experimental findings lfl3l . A similar negative MR was reported in 
Refs. irTTl [121 [T4ll . For ultra-high mobility samples [electron density 2 x 10 n cm~ 2 , vft~s ~ 80 pm, 
Tsm/Vs ~ 10, a ~ lOnm, ns ~ (2 pm)~ 2 \ 18]], one gets Bo ~ 60mT. A strong negative MR has indeed 
been observed |fT9ll29l in the very-high-mobility heterostructures, in qualitative agreement with the above 
theory. 
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2.3 Quasiclassical memory effects in ac magnetotransport 



In addition to the strong MR, the non-Markovian quasiclassical kinetics in the presence of long-range 
disorder gives rise to an anomalous ac response. In particular, the return-induced correction to the ac 
conductivity Recr(w) exhibits a kink ll30l oc \ui\ at u> — > 0. The quasiclassical zero-frequency anomaly is 
not sensitive to inelastic scattering (in contrast to the weak-localization quantum correction) but manifests 
itself only [ 3 1 ] in the presence of an external metallic gate that screens the long-range Coulomb interaction. 
As outlined below, in a strong B, the diagonal ac conductivity Kea xx (uj) shows pronounced resonant 
features [32 33] on top of the cyclotron resonance (CR), which are induced by the memory effects. 

In the Lorentz model (hard disks with no smooth disorder), which is intended to mimic a random 
AD array in heterostructures in the limit of a large spacer, the shape of the CR is very different from 
the Lorentzian and not characterized by the Drude scattering rate lf24l [321 . Altogether, the behavior of 
Re <t xx (lu) associated with quasiclassical cyclotron orbits skipping around the ADs turns out to be re- 
markably rich [32]. The skipping-orbit contribution is broadened on a scale of ui c and vanishes at uj c in a 
nonanalytical way as \u> — u> c \. Apart from these two features, Rc(t xx (uj) for moderately strong B with 
R c ^S> a oscillates with a period oj c up to u> = oj c R c /a and shows a series of square-root spikes for larger 
u>. The modulation yields exact zeros of the ac response at the harmonics of the CR. 

Adding a smooth random potential, present in typical heterostructures, changes the above picture in an 
essential way. Following the formalism of Sec. 12.21 the classical return-induced contribution to the real 
part of the oscillatory ac conductivity cr£f is represented for t$ <C T sm and 5 d as 11331 

4 C) /^ = -ReP^/nsrs , P w = £~ Jdte- iu p n [0,v F (t - 2nn/u c )] , (13) 



where a® is the real part of the dynamic Drude conductivity, Eq. (f3Qb - In Eq. (TT~3~b . the function p n (x± ,x\\) = 

(3a;^ + x 2 )/3n5 2 describes the electron distribution in time t = 2irn/ui c after n 



exp 



cyclotron revolutions along (2:11) and across (x±) the cyclotron orbit. In the case of a weak damping of 

(c) 

the quasiclassical magnetooscillations, cr^ is represented as a series of sharp dips at u> = Nlo c , whose 
amplitude and width are of order a D (a/NS)^^) 1 / 2 and t<7 1 + N 2 t^, respectively. In the regime of 



(c) 

harmonic oscillations (exponential damping), a u is written as [33| 

4 C) /^ = 1 - (a/7r 1/2 <5) cos(2ttlo/uj c ) exp [-(w/w c ) 2 (37r/w c T sm )] 



(14) 



It is worth stressing that these oscillations are of essentially classical origin and have nothing to do with 
the Landau quantization. The behavior of is illustrated in Fig. [5] where also the quantum oscillations 
oc exp(— 2ir/oj c T q ), with r q being the single-particle (quantum) relaxation time, are shown. One sees that 
the classical oscillations may be stronger than the quantum ones since in high-mobility structures r q <C T sm 
and the quantum oscillations are damped much more strongly. 



Q D 3 0.8 



0.4 
0.2 






(c) 





1/5 1/4 1/3 

CO /CO 



Fig. 5 Quasiclassical [oi , Eq. ( fT4t l and quantum [cri ] oscil- 
latory ac conductivity (normalized to the Drude conductivity <t°) 



vs ujc/lu for lo/2-k = 100 GHz, Tsm = 0.6 ns, 
T"s/r sm = 0.1, a/8 = 0.25 at uj c /lu = 1/2. 



i/Tq 



50. 
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3 Magnetotransport in modulated systems (lateral superlattices) 

3.1 Weiss oscillations in one-dimensional superlattices 

Transport properties of a two-dimensional electron gas (2DEG) subject to a periodic potential (lateral su- 
perlattice) with a period much shorter than the electron transport mean free path (but much larger than the 
Fermi wave length) have been intensively studied during the last decade. In a pioneering experiment [34| 
Weiss et al. discovered that a weak one-dimensional (ID) modulation with wave vector q || e x induces 
strong commensurability oscillations of the magnetoresistivity p xx (B) (while showing almost no effect on 
Pyy(B) and p xy (B)), with the minima satisfying the condition 2R c /a = n — 1/4, n = 1, 2, . . ., where 
R c is the cyclotron radius and a = 2tt jq the modulation wave length. The quasiclassical nature of these 
commensurability oscillations was demonstrated by Beenakker [35], who showed that the interplay of the 
cyclotron motion and the superlattice potential induces a drift of the guiding center along y axis, with 
an amplitude squared oscillating as cos 2 (qR c — 7r/4) (this is also reproduced by a quantum-mechanical 
calculation, see ll36l ). While describing nicely the period and the phase of the experimentally observed 
oscillations, the result of [35], however, failed to explain the observed rapid decay of the oscillation ampli- 
tude with decreasing magnetic field. The cause for this discrepancy was in the treatment of disorder: while 
Ref. [35 ] assumed isotropic impurity scattering, in experimentally relevant high-mobility semiconductor 
heterostructures the random potential is very smooth and induces predominantly small-angle scattering, 
with the total relaxation rate t^ 1 much exceeding the momentum relaxation rate r t ^ . The theory of 
commensurability oscillations in one-dimensional modulation, V(x) = r\Ep cosqx with 77 <C 1, in the 
situation of smooth disorder was worked out in [37). 

The starting point is the Boltzmann equation for the distribution function F(x, n) of electrons, 

CF{x, n) = — ew(a;)En ; C — v(x)nd r + tOcd^ — sin <fyu' {x)d f f > — C , (15) 

where n — (cos (/>, sin <fi) is the direction and v(x) — [2m(EF + eU(x))] 1 / 2 the magnitude of the Fermi 
velocity, and C is the collision integral. The resulting modulation-induced contribution to resistivity reads 
EO 



Apxx v2ql -Q-^-MQ)J-M- Q 



po 4 smhnp, ' qvFT q 



s 1/2- 

1 - ( 1 + ^Q 2 

Ttt 



(16) 



where we introduced the dimensionless parameter Q = qR c convenient to characterize the strength of the 
magnetic field. At low magnetic fields, Q Qdis> with Qdis = (29V 71 ") 1 ^ 3 ' tne oscillations are exponen- 
tially damped by disorder, and the magnetoresistivity saturates at the value Ap xx /p = r] 2 ql/A. (At still 
lower magnetic fields, Q > Q c h, with Q c h = 2/rj, and for a sufficiently strong modulation, rf/ 2 ql ^s> 1, 
an additional strong magnetoresistivity occurs, dominated by the channeled orbits, see Sec. 13.31 ) In strong 
fields, Q <C Qdis, the amplitude of oscillations increases as B 3 , 

&p xx /po = [(vqlf/vQ 3 ] cos 2 (Q - tt/4) . (17) 

In Fig. [6] the theoretical results are compared with experimental data of Weiss et al ll34l . The sample 
parameters are ll34l q = 27r/382nm, n e — 3.16 x 10 11 cm -2 , r tr = 52ps, the total relaxation rate is 
taken to be r^ 1 = (3ps) _1 . As is seen from the figure, with the modulation strength r) = 0.065 a 
very good description of the experimentally observed magnetoresistivity is obtained. (At lowest B, the 
experimental data show positive magnetoresistance discussed in Sec. 13.31 ) The difference between the 
long-range potential scattering and the isotropic scattering is illustrated in the right panel, where the results 
for the modulation-induced Ap xx are plotted for both models of disorder at the same value of r tl . As was 
shown in [38 1, modulation-induced commensurability oscillations can be also observed in attenuation and 
velocity change of a surface acoustic wave propagating near a 2DEG. 
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Fig. 6 Weiss oscillations in a ID superlattice. Left panel: Experimental data of Ref.|34| (dash-dotted line) compared 
with the theoretical results for the long-range potential scattering (full line). Parameters are: q = 27r/382nm, n e = 
3.16 x 10 11 cm -2 , t — 52ps, r s = 3ps, r\ = 0.065. Right panel: Grating-induced correction to the 2DEG reistivity, 
A.p xx /po, in units of tj 2 /4 for the isotropic potential scattering (r tr = 52ps, full line) and the long-range random 
potential scattering (rtr = 52ps, r q = 3ps, dash-dotted line). The sample parameters (q = 27r/382nm, n e = 
3.16 x lO^cm -2 ) are the same as in Ref.(54). 



In ||39l Weiss oscillations were studied for the vicintiy of the v = 1/2 filling of the lowest Landau level; 
see also related works ll40ll4D . Within the composite fermion theory, the problem is described in terms 
of fermions subject to a spatially modulated magnetic field and scattered by a random magnetic field. The 
magnetic character of modulation shifts the phase of Weiss oscillations, while the random magnetic fields 
considerably enhances their damping. The obtained results are in agreement with experimental studies 
ll42ll43ll . which confirms the validity of the composite-fermion description of the v = 1/2 state. 



3.2 Two-dimensional superlattices 

Magnetotransport in 2D superlattices with small-angle impurity scattering was studied in l44l . It was 
shown that the shape of the magnetoresistivity depends crucially on the parameter 7 = r) 2 ql/4. 

For small 7 (corresponding typically to a modulation strength not exceeding a few percent) the magne- 
toresistivity is given by the perturbative formulas ( MoT ) - ( fTTb (the same as in ID superlattices) up to the 
point Q ~ Qp = [O.lSfyql) 2 ] 1 ' 3 , where the correction Ap xx becomes of the order of the Drude resistivity 
po- For higher magnetic fields the Peclet number P ~ r\ql / (qR c ) 3 ^ 2 characterizing the advection-diffusion 
problem becomes large and the transport is determined by a narrow boundary layer around a square net- 
work of separatrices. As a result, the B 3 -dependence of the oscillation amplitude characteristic for the 
perturbative (Q > Qp) regime crosses over to a much slower _B 3 / 4 -increase at Q < Qp, 

Ap xx /p = (8ny/ 4 C( m l)^ 2 Q- s ^\ cos(Q - ir/i)^ 2 . (18) 

For 7 ^> 1 (which is typically valid for the modulation strength 7/ larger than 10 -j- 15%) the oscillations 
are damped at low magnetic fields not by disorder (as in the perturbative regime) but by the modulation- 
induced chaotic diffusion. The oscillations become observable at Q ~ Q a d = (47r /77 2 ) 1 / 3 where the 
motion of electrons in the superlattice potential acquires the form of adiabatic drift. Since the violation of 
adiabaticity is exponentially small, the magnetoresistivity drops exponentially in a logarithmically narrow 
interval of magnetic fields, Q' ad = Q a d(ln7)~ 2 ^ 3 < Q < Q a d, 

Ap xx /p cx exp[-(^/2V2| cos(Q - 7r/4)|)(Q ad /Q) 3 / 2 ] . (19) 
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At higher magnetic fields the impurity scattering starts to dominate over the non-adiabatic processes and 
thus to determine the diffusion constant of the advection-diffusion problem, so that the commensurability 
oscillations take the same form ([T8l as in the large-P limit of the 7 <C 1 regime. 




! di S b p b Bad B = d B 



Fig. 7 Left panel: Schematic representation of the magnetoresistivity Ap xx (B) induced by a 2D modulation in the 
case 7< 1. Characteristic points Bdis and Bp on the magnetic field axis (corresponding to Q = Qdis and Q = Qp, 
see the text) are shown. Below £>di s the oscillations are exponentially damped and the magnetoresistivity saturates 
at Ap xx = 7, while at B = Bp the _B 3 -behavior of Ap xx changes to a much slower, B 3 / 4 -increase. Right panel: 
Schematic representation of Ap xx (B) in the case 7^1. The magnetoresistivity starts to drop exponentially and the 
commensurability oscillation appear at the value B a( j of the magnetic field where the motion in the periodic potential 
takes the form of an adiabatic drift. At B ~ B ad the disorder starts to dominate over the non-adiabatic effects, leading 
to a _B 3//4 -increase of the oscillation amplitude. 



3.3 Low-field magnetoresistance 

A distinct low-field magnetoresistivity was observed, along with the commensurability oscillations, in the 
original experiment |34|, as well as in numerous later experiments on the transport in a lateral superlattice. 
Specifically, in low magnetic fields B a positive magnetoresistivity was found, followed by a maximum 
in Pxx(B). For not too strong modulation, the relevant magnetic fields are much weaker than those where 
the Weiss oscillations are observed, so that the two effects can be easily separated. Soon after the first 
experimental observation it was understood [47 1 (see also [46 1) that the low-field magnetoresistivity is 
related to the existence of open (channeled) orbits in the magnetic fields B < B c = (rjc/2e)qmVF. It is 
worth mentioning that this effect, which is not found within the 77-expansion used in Refs. [35 , 37], has its 
counterpart in the context of the sound absorption in metals in the presence of a magnetic field. There, the 
trapping of electrons in channeled orbits by a sound wave leads to non-linearity of the acoustic response of 
an electron gas, as was observed experimentally [48 1 and analyzed theoretically ||49l . 

A quantitative analytical description of the problem in the presence of disorder was worked out in [45 1. 
It was found that for a suffuciently strong modulation, rf/ 2 ql 1, the contribution of channeled orbits to 
resistivity has the form 

Apfjpo = (V2/7T 2 W /2 (ql) 2 Fc h (P) , (20) 

where (3 = B/B c and -F c h(/3) is a parameterless function shown in the left panel of Fig. [8] This induces a 
low-field magnetoresistance that scales as rf/ 2 with the modulation strength. It was further shown in [45 1 
that the contribution of non-channeled orbits is also modified at B < B c , 

Ap™/p - 2klu c t(v 2 )/v 2 f = ( V 2 /27r)qlF ac (B/B c ) , (21) 

where the dimensionless function F nc ((3) is shown in the middle panel ofFig.[8] In the right panel ofFig.[8] 
the theoretical results are compared with experimental data of Ref. 11501 . 
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B P B (T) 



Fig. 8 Left panel: Function F c h({3) describing magnetic field dependence of the contribution of channeled orbits 
to the resistivity. Middle panel: Function F nc (f3) characterizing the magnetic field dependence of the contribution of 
non-channeled orbits to resistivity. The dashed line indicates the asymptotic value F nc (f3 2> 1) = n/2. Right panel: 
Experimental data for the low-field magnetoresistivity from Ref. [50| (dashed curve) compared to the theoretical results 
for the contributions of channeled rEq. (I21l>1 and non-channeled rEq.d20H orbits. The characteristic field B c and the 
modulation amplitude found from the fit are B c ~ 0.37 T and r\ ~ 0.16. 



4 Interaction-induced magnetoresistance in non-quantizing fields 

As discussed in Sec. [2] the longitudinal resistivity of an isotropic degenerate system is £?-independent 
within the Drude-Boltzmann theory, p xx [B) = po = (e 2 z/uf,7t r )~ 1 , where Vq is the density of states per 
spin. There are several distinct sources of a non-trivial MR, which reflect the rich physics of 2D systems. 
First, quasiclassical memory effects may lead to a MR(see Sec|2]), which shows no T-dependence at low 
temperatures. Second, weak localization BTI induces a negative quantum MR restricted to the range of 
weak magnetic fields. Another quantum-mechanical source of the MR is electron-electron interaction. 
This type of the MR is the subject of the present Section. 

It was discovered by Altshuler and Aronov IfSTl that the Coulomb interaction enhanced by the diffusive 
motion of electrons gives rise to a quantum correction to conductivity, which has in 2D the form 8a xx ~ 
(e 2 /27r 2 ) In Tr tt (it is assumed here for simplicity that n -C kp, where n — Aire 2 ^ is the inverse screening 
length; we also set ks = fr = 1). The condition Tr tl -C 1 under which this result is derived BTI implies 
that electrons move diffusively on the time scale 1/T and is termed the "diffusive regime". Subsequent 
work [52 1 showed that in a strong magnetic field this correction (in combination with Sa xy = 0) gives rise 
to a parabolic interaction-induced quantum MR, 

S Pxx {B)/p ~ [(w c r tr ) 2 - IKttM)" 1 lnTrtr, Tr tr « 1, (22) 

where lo c — eB/mc is the cyclotron frequency and / = vpr tl - the transport mean free path. 

The effect of interaction on the conductivity in the "ballistic regime" T > 1 /r tl has attracted a great deal 
of interest in a context of 2D systems showing a seemingly metallic behavior, dp/dT > [53, 54] . Zala, 
Narozhny, and Aleiner lf55ll developed a systematic theory of the interaction corrections valid for arbitrary 
Ttm. In the ballistic range of temperatures, this theory predicts a linear-in-T correction to conductivity a xx 
and a 1/T correction to the Hall coefficient p xy /B at B — > 0, and describes the MR in a parallel field. 
The consideration of [55] is restricted, however, to classically weak transverse fields, uj c T a <C 1, and to the 
white-noise disorder. 

In this Section, we present a general theory of the interaction-induced corrections to the conductivity 
tensor of 2D electrons valid for arbitrary T, B and type of disorder ||56ll . A general expression for 8a a p 
is derived in terms of the ballistic propagator D(u>, q; n, n') describing the quasiclassical propagation of 
an electron in the phase space (n is the unit vector characterizing the direction of velocity on the Fermi 
surface). The result for the exchange contribution reads 

Sa af3 = -2e 2 vlv J |w coth^j ^ Im [ (7(w, q) q) ] , (23) 
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where U(ui, q) is the interaction potential equal to a constant Uq for point-like interaction and to U(ui, q) = 
2iTe 2 {q + k[1 + iuj(D(ui, q))]}" 1 for screened Coulomb interaction. The angular brackets (...) denote 
averaging over velocity directions n, n'. The tensor B a p(u>, q) is given by 



B a p(u,q) = T cer Tru [(DS J eD)-2(Dn J Wn e D)]T B i3 
— 2T ai {n~ l DnpD) — (Dn a DnpD), 



T ai (5 jS {D)/2- (n 7 Dn s ))T S p 

(24) 



where T afj = 2 {n a Dnp)\ q=Q ^ = <x a p/e 2 v F v , S xx = S yy = W(n,n'), ^ xy - -^ yx 
u c /2wvq, and W(n, n') is the impurity scattering cross-section. At B — > one recovers the results 
for ^(Tica; and obtained in a different way in [55 1 for a white-noise disorder. Needless to say, in the 
diffusive limit, Eqs. d23T >. (l24l i reproduce (for arbitrary B and disorder range) the logarithmic correction 
(ES. 

The structure of Eqs. (l23l l. (124-b implies that the interaction correction is governed by returns of a particle 
to the original point in a time t < T^ 1 <C r a . In a smooth random potential with a correlation length d 
fc^ 1 the return probability is exponentially suppressed for t <C r tr . Therefore, the interaction correction in 
the ballistic regime is exponentially small at B — for the case of smooth disorder. Moreover, the same 
argument applies to the case of a non-zero B, as long as uj c -C T. 

The situation changes qualitatively in a strong B, lu c ^> T, t^ 1 : the particle experiences within the time 
t ~ T _1 multiple cyclotron returns to the region close to the starting point. The MR is then determined by 
the correction to a xx . For the Coulomb interaction, the exchange contribution to the MR is given by 



5Pxx/P0 = -{UcT*) 2 G F {TT ti )llTk F l. 



(25) 



For the point-like interaction a similar result is obtained, with the replacement Gf{Ttx x ) — * i^oC/oGo(7 1 T tr ). 
The functions Go(Tr tr ) and G , F(7 1 T tr ), governing the T-dependence of the MR, are shown in Fig. [9] In the 
diffusive (Tr tl - -C 1) and ballistic (Tr tr ^> 1) limits they have the following asymptotics 



G (x) 



In a; + const, 



c x 



-1/2 



where c = 3C(3/2)/16VtF ~ 0.276 





G„(x) 




1. 7-111 X 




0.276 X"" 2 






a \ 





X < 1, 

x > 1, 





G F (x) 




1.6-ln x 




0.138 x~ ,a 


b 





T- 



Fig. 9 Functions G (Tr tr ) (a) and G F (Tr, r ) (b) deter- 
mining the T-dependence of the exchange term for point- 
like and Coulomb interaction, respectively, Eq. ( I26l >. 



G F (x) 



— Inx + const, 
co*- 1/2 A 



a; < 1, 
x~> 1, 



(26) 




Fig. 10 Functions Go lix (Tr tr ) (a) and G^ ix (Tr tr ) (b) de- 
scribing the T dependence of the MR for point-like and 
Coulomb interaction, respectively, in the mixed-disorder 
model for different values of parameter 7 = r sm /r tr = 
20, 10, 5 (from top to bottom). Dashed curves represent 
these functions for purely smooth disorder (7 = 1). 



We turn now to the Hartree term, assuming n -C hp- The expression for its triplet part is analogous to 
d23l l with the replacement of U (to, q) by — |J7(0, 2k F sin[(0 — </>')/2]), where <fi and <f)' are the angles of 
the electron velocity. As to the singlet part, it is renormalized by mixing with the exchange term. The total 
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Hartree contribution (see Fig. 2) reads 



r.H /m I V In 2/ [f ln(Tr tr ) + lny] , Tt u < 1, 

PxxK } = frac{u c T«)\ 2 k F l I yln 2 ^^-) 1 / 2 ], 1 « Tr tr « (fc F / K ) 2 , (27) 



= frac{uj c r a f^k F l { y\n 2 [y{l 
P0 { nco(TT*)- 1 / 2 , Tr tr » (fc F / K ) 2 . 



If K/fcp is not small, the exchange contribution remains unchanged, while the Hartree term is subject to 
strong Fermi-liquid renormalization 0111551 and is determined by angular harmonics F°f of the Fermi- 
liquid interaction F a ' p (9). The theoretically predicted T -1 / 2 dependence of the interaction-induced MR 
has been observed in a high-mobility n-GaAs heterostructure with the smooth disorder 11571 . 

Calculation of the correction 5p xy to the Hall resistivity requires evaluation of both 8a xx and 5a xy . The 
temperature dependence of 8p xy in a strong B is governed by 5a xx in the diffusive limit and by 8a xy in 
the ballistic limit. For the point-like interaction, we find 

5p xv /p xv = p U Q G p ^(Tr tx )/7rk F l, G%"(z) = l l]^,^ ! f ]' (28) 



-llCi2; 1/2 , Z>1, 

with c\ = — y / 7rC(l/2)/4 ~ 0.647. An analogous consideration for the Coulomb interaction yields a 
similar result for the exchange correction (Fig, [nil 

5pl y /p xv = Gp v W/irW, G£" (x) = { I^^V?' Itl. (29) 



CD 



5 

-5 
-10 
-15 



-20 ' — 
0.01 



-3.56 x -0.9 



0.1 1 

Tx 



Fig. 11 Functions Gq WV (Tt 1t ) (lower curve)and Gp (Tr a ) (upper curve) de- 
scribing the temperature dependence of the Hall resistivity for point-like and 
10 Coulomb interaction, respectively. Diffusive (a; -C 1) and ballistic (x 2> 1) 
asymptotics are also shown. 



Above the interaction correction for a system with a small-angle scattering induced by smooth disorder 
with correlation length d ^> kp 1 has been studied. This is a typical situation for high-mobility GaAs 
structures with sufficiently large spacer d. It is known, however, that with further increasing width of 
the spacer the large-angle scattering on residual impurities and interface roughness becomes important 
and limits the mobility (see Sec. 12.21 ). Furthermore, in Si-based structures the transport relaxation rate is 
usually governed by scattering on short-range impurities. This suggests considering the two-component 
model of disorder: white-noise random potential with a mean free time r wn and a smooth random potential 
with a transport relaxation time r sm . It is assumed that while the transport relaxation rate t,^ 1 = t",^ + 
Tg^ 1 is governed by short-range disorder, r wn <C r sm , the damping of SdHO is dominated by smooth 
random potential. This allows one to consider the range of classically strong magnetic fields, w c r wn ^> 1 , 
neglecting at the same time Landau quantization. 

The presence of short-range scatterers enhances the MR as compared to the case of smooth disorder. 
For the point-like interaction the MR reads 



5p xx (B) _ (LU c T tI ) 2 V U 



Po 



irkpl 



Gmix 




rri T Sm 

T T tt , 



Gmix 




— In a; -I- 
4c o7 1/2 . 



(27) 

„-l/2 



1/2 



X < 1, 
X > 1. 
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In the case of Coulomb interaction, the exchange contribution is given by 

6o F > mi *(B) (cut,) 2 ( r \ { - In a; + h/2) 1/2 , x < 1. 

"/'.,., U> , Mr I ^mix I rp T T sm \ Q mix ( x 7) ~ ' 



Po irk pi V Ttr / 1 2c 7 1/2 x" 1 / 2 , x>l. 



These results (illustrated in Fig. 1 1 Ob are in a good agreement with experimental data of Ref. [58 1 for mixed 
disorder in a Si/SiGe heterostructure. Finally, the ballistic contribution to p xy is also enhanced by a factor 
~ (TsmAtr) 1 ' 2 in the mixed disorder model. 

The formalism can be further applied to anisotropic systems. The interaction-induced correction mixes 
the components p xx and p yy of the resistivity tensor. This result is of special interest in the case of sys- 
tems subject to a one-dimensional periodic modulation (lateral superlattice, Sec. [31 wave vector k || e x ). 
Specifically, it has been shown that the interaction induces novel oscillations in p yy , which are in phase 
with quasiclassical commensurability (Weiss) oscillations in p xx . 



5 Influence of Landau quantization on magnetotransport 

In this section, we address magnetooscillations in the dissipative dc and ac conductivity of a 2DEG gov- 
erned by the Landau quantization. Despite these effects, to the first place, Shubnikov-de Haas oscillations 
(SdHO), are well established experimentally, the theoretical description until recently was only available 
[59 1 for fully separated Landau levels (LLs) with point-like scatterers ll60l . A systematic approach to the 
problem was developed in Ref. ll6D . The results of this work, valid also for overlapping LLs, and for 
experimentally relevant case of smooth disorder, with the correlation length d^$> kp , are reviewed below. 

Within the quasiclassical Boltzmann theory, the dissipative ac conductivity <j u = a + (oj) + <J- (lu) of a 
non-interacting 2DEG is given by the Drude formula (we neglect spin for simplicity), 

o£(w) = (l/4)eVo4V[l + K ± ufrl] , (30) 

where Vq = to/27t and r tr are the density of states (DOS) and the transport relaxation time at B = 0, 
uj c = eB/mc the cyclotron frequency, and m is the electron effective mass. We consider a 2DEG sub- 
jected to quantizing magnetic field B and a random potential U (r) characterized by a correlation function 
(U (r)U(r')) — W(\r — r'|). The total and the transport relaxation rates at B = are 

\\ } = 2-kvo I ^iy(2fc F sin^) x ( I , 
T t / / °J 2ir v 2 J \ (1-COS0) 

where VK(q) is the Fourier transform of W(r). While we are mainly interested in the experimentally 
relevant case of smooth disorder, d ^> kp 1 , with T tr /r q ~ (kpd) 2 S> 1, our results are valid for arbitrary 
d (i.e., including short-range disorder with T tr /r q ~ 1). The conductivity is given by the Kubo formula 

a u = - {e 2 /^Vuj)jde (f E - f £+ul )Tr v x (G^ +LU - G? +LU )v x (G? - Gf) , (31) 

where f e is the Fermi distribution, G R > A are the retarded and advanced Green functions, the bar denotes 
impurity averaging, and V is the system area. At high LLs, ep ^S> ui, uj c , disorder can be treated within the 
self-consistent Born approximation (SCBA) [59| provided the disorder correlation length satisfies d <C Ib 
and d <C vpT q , where Ib = (c/eB) 1 / 2 is the magnetic length (62]. The SCBA equations for the Green 
function in the LL representation, G R — (G£)*, read 15911621 , 

G%{e) = (£-£„- S,)- 1 , S e = (w c /27rr,, ) G%(e) , (32) 

where e n = (n + ^)u c is the n-th LL energy (Fig. [T2k). The conductivity OTb is given by an electronic 
bubble with a vertex correction, i.e., by a sum of ladder diagrams, Fig.fTZb.c. In the case of white-noise 
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disorder, r q = T tr , the vertex correction is zero, and it suffices to evaluate the bare bubble 



4M= 



e 2 VpV f de 



J ^ 



Re (IT 



RA_ U R R) 



n 



RR(RA) _ U c ^ 



± 



For the case of smooth disorder we have to take into account the vertex correction (Fig. lc) while averaging 
in Eq. OTb . This is a non-trivial task since the disorder mixes strongly the LLs, thus seriously complicating 
a direct calculation in the LL representation. The result, however, acquires a remarkably simple and phys- 
ically transparent form: provided the above SCBA conditions are fullfilled, the inclusion of the vertex cor- 



rection results in a replacement of jj^^" 4 ^ above by Tl 1 ^. 11 ^ = 

It follows that has a Drude-type structure with the DOS, v{e) 
transport time renormalized due to Landau quantization: 



n 



RR(RA)\ { i : 
7r _1 ^ow c Im^ n G„(e), and the 



a± (uj) 



2„,2 



e v 



de (f e - fe+u,) v{e) T tr p{e + u) 
[T u 2 B (e) + T-%(e + w)]/2 + (w ± u c ) 



Ttr Me) 



Ttr 



v{e) 



(33) 



Formula d33l is the main result of this section. Let us emphasizethat the single-particle time r q enters 
Eq. (|33l l only through the DOS; everywhere else it has been replaced by the transport time r tr due to the 
vertex correction. In the following, we analyze Eq. d33l in several important limiting cases. 



(a) 





Fig. 12 (a) SCBA equation for the Green function; (b) 
dynamical conductivity with vertex correction (c). 



(0/<B 



Fig. 13 Magnetooscillations of the dynamical 
conductivity for r tr /T q = 10. Solid line: sepa- 
rated LLs, uj c T q /-K = 3.25; dashed line: over- 
lapping LLs, uj c T q /ir = 1. Inset: a xx for fixed 
uiTq/2-K — 1 as a function of u c . 



In the regime of strongly overlapping LLs, ui c T q 1, the solution to SCBA equations d32b is most 
easily obtained using the Poisson formula, — Ylk I dxF(x) exp(2irikx). The k = term yields 

the B = result, while the k = ±1 contributions provide the leading oscillatory correction to the DOS, 



v{e) = v [l - 25cos(2ire/bJ c ) + 0{8 2 )] , 8 = exp(-7r/w c r q ) < 1 
To first order in 8, Eq. d33l produces the following result: 



n , z =1-28 T \ I cos 

cr^(cj) 



2a \ sin(27T£j/a; c ) 3qj_ + 1 sin 2 (7rtj/ijj c ) 



+ 1 2t:lu/uj c 



+ 1 a±irLu/uj c 



(34) 



, (35) 



where a± = r tr (w ± oj c ), and the Dingle factor J-(X) = X/ s'mhX describes the T-damping of the 
SdHO. In the dc limit to — > 0, this result confirms the form of SdHO in smooth disorder conjectured in 
[63 1. If T is higher than the Dingle temperature To = l/2nT q , the temperature smearing becomes the 
dominant damping factor. In high-mobility 2DEG the Dingle temperature is as low as Tjj ~ 100 mK, so 
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that for characteristic measurement temperatures T ~ 1 K the first-order correction d35l l will be completely 
suppressed. However, there exists a correction of order S 2 , oscillatory in u>/ lu c , which is not affected by the 
temperature. To obtain it, there is no need to calculate v(e) to second order, since the corresponding terms 
oscillate with e, and doesn't survive the high-T limit. The leading quantum correction at T ^> To results 
from the averaging Eq. (1331 over fast energy oscillations of the first-order v {£), Eq. d34l >. which gives 



a£ 2) H = a£(w)<j 1 + 2(5 



a 2 ±(a 2 



3) 2lTUJ 
77T- COS 



0J C 



a±(3ai - 1) 2vrcj 
~TyT~ sm 



(36) 



The regime which is most interesting theoretically and relevant experimentally is that of long-range 
disorder, r t i /r q ^> 1, and a classically strong magnetic field, uu Cl w ^> r t 7 . In this situation Eq. d3~3b reads 



j Is fe+u; i \ / \ 

de 2 — u \ £ ) v \ e + w ) j 



E 



e 2 is vj 



4T tr (w ± uj c ) 2 



(37) 



or, in dc limit, 



-agl de[v 2 {e)/v 2 ]d e f(e) , a£ = e 2 ^ oU |/2r tr ^ c 2 , (38) 
In the limit of separated LLs, uj c r q 1, the DOS is a sequence of semicircles of width 21" <C u c , 

(39) 



v(e) = ^r q ^ n ReVr 2 - (e - e„) 2 , T - ^2lu c /tt T(1 . 

In this case, <r w is non-zero only for uj in intervals [Afw c — 2r, Muj c + 2T] with an integer M. 

Oscillations in a u with lo/lo c for 2DEG with smooth disorder at w c r tr 3> 1 and T ^> To are illustrated 
in Fig. [13] For overlapping LLs, oscillations away from the cyclotron peak are described by simple formula 
(7^/(7° = 1 + 28 2 cos(27rw/o; c ). For separated LLs, at the center of the M = 1 interval we find a CR 
peak of height a xx (u> = u) c ) = (e 2 i'ov 2 7 /7rr)r tri o/T q and width ~ rr q /T tr ,o. All other peaks (M ^ 1) are 
smaller by a factor ~ T tr /r q ~ (kpd) 2 3> 1, 



£7 ra (w = Mlo c ) = (4e 2 t,o4r/37r W 2 )(r q /r tr )[(M 2 + 1)/(A/ 2 - l) 2 



(40) 



6 Interaction effects in quantizing magnetic fields 

6. 1 Interaction effects on oscillations 

In this Section we discuss the interaction effects on magnetooscillations (de Haas-van Alphen and Shub- 
nikov - de Haas oscillations), closely following Ref . [ 64 1 . In addition to experimental motivation related to 
the apparent metal-insulator transition in two-dimensional systems, the corresponding theory complements 
the recently developed theory of interaction effects in transport of 2D electrons in zero and non-quantizing 
magnetic fields Il55ll56l . 

The starting point is the expression for the thermodynamic potential derived in the paper by Luttinger 
and Ward (55) 

n = -TTrln(-G _1 ) - TTr(GS) + fi', (41) 

where the trace implies summation over Landau levels N and over fermionic Matsubara frequencies e n = 
[2n + l)wrT, G(ie n ,Nuj c ) = [ie n + fi — (N + l/2)w c - E(ie„, iVw c )] _1 is the dressed Matsubara 
Green's function, and S(ie n , Nui c ) is a self-energy part of Green's function which includes all the disorder 
and interaction effects. The terms — TTr(GS) and fi' in fiTt are introduced to avoid double-counting of 
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(a) 



y 



Fig. 14 Self energy diagrams in the first order in the effective interaction 
(wavy line). Black triangles denote impurity ladders F dressing interac- 
tion vertices, dashed line is a single-impurity line. Diagrams (b) repre- 
sent the "Hikami-box" contribution to the self-energy which restores the 
gauge-invariance of the damping. 



diagrams [65, 66 1. The term SI' denotes the sum of all so-called skeleton diagrams with all bare Green's 
functions replaced by dressed Green's functions. 

As shown in Ref. 11671 . the exponential decay of magnetooscillations is described by the Trln-term. 
The oscillatory parts of the additional terms, which are introduced to correct for overcounting, cancel each 
other. In order to obtain the correction to the thermodynamic potential one needs to calculate the self- 
energy part of the Green's function. It is worth noting that the inelastic contribution to the self-energy oc 
[(ttT) 2 — e^]signe„ vanishes for e n — ttT, and thus does not affect the damping of the magnetooscillations 
B(T) for T > lu c , in agreement with Refs. Ii68ll69l 

For overlapping Landau levels [Eq. d34lil. the magnetooscillations are damped by disorder even at zero 
temperature via the standard Dingle factor 6 — exp(— 7t/w c t 9 ). Therefore we will consider only the 
first harmonics of the oscillations, A\, neglecting all higher harmonics whose damping is much stronger. 
In what follows we concentrate on the case T lu c . Under this condition, the first harmonics of the 
oscillatory part of the thermodynamic potential 

il osc ~ 2v (lo c /2tt) 2 A 1 cos(2ir 2 n e /eB), (42) 

is given by 

A 1 = (4tt 2 TK) exp [~2n 2 T/Lu* c - tt/uj*^] exp[B(T)], (43) 
which is a standard FL Lifshitz-Kosevich expression multiplied by the additional factor with 

B(T) = -2TriZ5T,(iTrT,£ a )/uj*. (44) 

Here lo* = eB / m* is the Fermi-liquid (FL) renormalized effective cyclotron frequency in a pure system 
at zero T, which is related to the FL-renormalized effective mass m*, Z is the FL Z-factor (given by 
the residue of the Green's function), r* is the FL-renormalized scattering time, and 5£(ie n ,£o) is the 
self-energy part (taken at the pole £o of the Green's function in the presence of disorder) describing the 
interplay of disorder and interaction. 

Hereafter white-noise disorder with r tl = r q = t is considered. Evaluating the sum of six digrams 
(Fig. H~4b for <5E, one gets the following expression for the damping exponent in the case of short-range 
interaction £/ : 

B{T) - -const v U (ir/Lu c T q ) + (ttT/uj c ) {v U /s f t) ln(e F /T). (45) 

The first term in Eq. d45l ) describes the T-independent FL-renormalization of r q due to vertex corrections 
and should be included in the effective relaxation time r* . The second term represents the T-dependent 
contribution to the damping factor that we are interested in and is analyzed below. 

The above result d45l > can be interpreted in terms of corrections to the effective mass (or u> c ) and the 
quantum elastic scattering rate r q entering the standard Lifshitz-Kosevich formula. These corrections come 
from the interplay of disorder and interaction, leading to 

B(T) = -(2n 2 T/uj c )(Sm/m) - (tt/lu c t q ) [6m/ m - 5 Tq /T q ] . (46) 

It is worth noting that the FL-renormalization does not affect the product ui c m = eB. 
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Comparing d45b and ( f46b , one can see that the T In T dependence of the damping factor could in prin- 
ciple originate either from the In T correction to the effective mass, or from the T In T-type correction to 
T q . This led the authors of Ref. |68| to the conclusion that the nonlinear T— dependence of the damping 
factor may be equivalently interpreted either as a T— dependent renormalization of the effective mass or 
as a T— dependent Dingle temperature. It is clear, however, that these two possibilities correspond to dif- 
ferent physical processes. To identify the physical origin of the leading contribution to the damping it is 
instructive to obtain B(T) using the expression for the self-energy analytically continued to real values of 
energies e n — > — ie . 

Having calculated ReS and ImE for real energies e, one can determine 6m and Sr q . Indeed, the magni- 
tude of the first harmonics of the magnetooscillations of the thermodynamic density of states is expressed 
through the real-e self-energy <5S(e) as follows: 



A 1 (T) = -J deAiierfdeMe), 



At (e, T) = exp { — [e - Re<5£(e, £„)] } cxp ( — + — Im5E( £ , &) 

L0 r " UJrTn 0J r 



exp 



2nie 



Sm(s, T) 



exp 



Jc'q 



1 



c ' q 

Sm(e,T) S T q(s,T) 



(47) 
(48) 
(49) 



where /t(e) = [1 + exp(e/T)] 1 is the Fermi distribution function. This allows one to express <5r 9 (e, T) 
and 5m(e, T) through ReE(e) and ImS(e) as follows: 



Sm(e, T)/m 

STq(e,T)/T q 



-e _1 Re<JE(e, T) = -{v U Q /2%e F T)\n{E F /max[\e\,T}}, 
2T q ImSY,(e,T) + 5m(e,T)/m 



T 

vqU — In 
e F 



In 



e F 



2ire F T max[|e|,T] 



(50) 



(51) 



It is clear from these results that the leading term in B(T) [proportional to T\n(e F /T), Eq. (05])] origi- 
nates from the real part of the self-energy, i.e. from renormalization of the effective mass, which affects 
incommensurability of the oscillations at different values of energy e. The contribution of the imaginary 
part of the self-energy, which is governed in the ballistic regime by the renormalization of the scattering 
time, is smaller by a factor \n(e F r). The obtained result for the interaction-induced correction to the quan- 
tum scattering time r q , Eq. (IBTT l. agrees, up to a factor |, with the correction to the transport time following 
from the calculation of conductivity correction in the ballistic regime in Ref. If55ll . 

In the case of Coulomb interaction, one should take into account the dynamical screening of the interac- 
tion within the random phase approximation (RPA). This leads to different asymptotics of the self-energy 
in the diffusive and ballistic regimes, in contrast to the case of weak short-range interaction. The In- 
dependence of the leading correction to the magnetooscillations damping factor due to the interaction in 
the singlet channel has the form (Fig. IT31 



B P (T) 



7T T 

LU c T q £p 



(3/2)ln(e F /T) - (1/2) ln(47re F r), AttTt < 1, 



(52) 



ln(e F /T), 



4itTt > 1. 



Calculation of the corresponding triplet contribution leads to qualitatively similar asymptotics. The leading 
term in the total correction to the damping factor in the ballistic regime, realized in experiments on low- 
disorder samples at realistic temperatures, takes the simple form 



B(T) = B P (T) + B°(T) ~ [1 + 3F CT (1 + i^)" 1 ] (nT /to c T q e F ) \n{e F /T). 
As discussed above, this result arises due to the correction to the effective mass. 



(53) 



© 2003 WILEY- VCH Verlag GmbH & Co. KGaA, Weinhcim 



18 



LA. Dmitriev et al.: Magnetotransport in quantum-Hall systems 



3 



0.5 




a 



b 



-1 



-0.5 







2 







0.1 



Tt 



Tx 



Fig. 15 Temperature dependence of the singlet channel correction to the damping factor B P (T) for Aire ft = 100 
(solid line) with the low-T (dot-dashed) and high-T (dashed) asymptotics. (a) Wide temperature range: on this scale 
B P (T) is essentially indistinguishable from its high-T asymptotics; (b) low-T part: the crossover between the two 
asymptotics occurs at Tt ~ 0.05. 

6.2 Coulomb drag in high Landau levels 

Coulomb drag between parallel two-dimensional electron systems IfTOlrTTl has developed into a powerful 
probe of quantum-Hall systems f72l |73l |74l |75l |76l 1771 [78l [79l , providing information which is com- 
plementary to conventional transport measurements. The drag signal is the voltage V developing in the 
open-circuit passive layer when a current I is applied in the active layer. The drag resistance (also known 
as transresistance) is then defined by Rr> = V/I. As a function of interlayer spacing a, the interlayer 
coupling changes from weak at large spacings where it can be treated in perturbation theory, to strong at 
small spacings where it can result in states with strong interlayer correlations 17611771 . 

In a simple picture of Coulomb drag, the carriers of the active layer transfer momentum to the carriers 
of the passive layer by interlayer electron-electron scattering. The phase space for interlayer scattering is 
proportional to the temperature T in either layer predicting a monotonous temperature dependence Rr> oc 
T 2 of the drag resistance. Moreover, the signs of the voltages in active and passive layer are expected to 
be opposite (the same) for carriers of equal (opposite) charge in the two layers [80). 

Remarkably, experiments show that Coulomb drag behaves very differently from these simple expecta- 
tions when a perpendicular magnetic field B is applied such that the Fermi energy £p is in a high Landau 
level, £p/hjj c ^> 1. (oj c is the cyclotron frequency.) Several experiments 117411781 in the regime of weak 
interlayer coupling observed negative drag when the filling factors in the two layers are different. A more 
recent experiment [79] also reveals a non-monotonic dependence on temperature. While the drag resis- 
tivity shows a quadratic temperature dependence at sufficiently high temperatures, where drag is always 
positive, an additional peak develops at low temperatures which can have both a positive or a negative sign 
depending on the filling-factor difference between the two layers. 

In this Section, we present the theory of Coulomb drag in the limit of high Landau levels l82l . In 
a strong magnetic field, u> c T tl - ^ 1, the intralayer Hall resistivity p xy dominates over the longitudinal 
resistivity p xx . Therefore, the drag resistivity is given by 



The Coulomb drag in strong magnetic fields is an interplay of two contributions, as illustrated in Fig. [TBI 
At high temperatures, the leading contribution is due to breaking of particle-hole symmetry by the curvature 
of the zero-B electron spectrum. This "normal" contribution to the drag is always positive and increases 
in a broad temperature range as T 2 . At low temperatures, another, "anomalous", contribution dominates, 
which arises from the breaking of particle-hole symmetry by the energy dependence of the density of 
states related to Landau quantization. This contribution is sharply peaked at a temperatute T ~ A (where 



Pxx — 



Pxy r yy Pyx 



(54) 



© 2003 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



pss header will be provided by the publisher 



19 



p-h symmetry p-h asymmetry 

curvature LL - quantization 



p-h symmetry 



E / 








k 




ph i/h 






v = dE/dk 







p-h asymmetry 

DoS V(E) 



•- 



o- 

h O- 



O- 



o 
o 
o- 
o 



J p. h 



Fig. 16 Schematic illustration of different sources of particle-hole asymmetry: curvature of zeio-B spectrum E(k) 
vs LL-quantization of the density of states (DoS) u(E). In the particle-hole (p-h) symmetric case, the electronic and 
hole contributions to the current induced in the passive layer (j e and jh , respectively) compensate each other. When 
the p-h asymmetry is generated by a finite curvature, the velocities of electrons and holes (shown by arrows in the right 
panel) are different, which destroys the compensation. This is the "conventional" mechanism of the drag. When the 
DoS depends on energy (in the present case because of the LL-quantization), an "anomalous" drag arises due to the 
difference in numbers of occupied electronic and hole states. 




A = 2r is the Landau level width) and has an oscillatory sign depending on the density mismatch between 
the two layers. 

Since the momenta transferred from one layer to the other are effectively restricted by the inverse in- 
terlayer distance, a , the behavior of the transresistivity will be essentially dependent on the relation 
between R c and a. Specifically, with increasing R c /a the following four regimes are identified i) diffusive, 
R c /a < 1, ii) weakly ballistic, 1 < R c /a < to c /A, iii) ballistic, uj c /A < R c /a < NA/uj c , and iv) 
ultra-ballistic, NA/uj c <§; R c /a. In all regimes, the temperature-dependence of the drag resistivity is non- 
monotonous: the absolute value of p^ x (T) shows a peak around T ~ A and increases again at T ^ uj c . 
However, the T— and B— dependences of p® x , as well as the sign of the low-temperature peak (the high- 
temperature drag is always positive), are specific for each particular regime, as illustrated in Fig. [771 and 
summarized below. 

Diffusive regime, R c /a <C 1. In the diffusive regime, the drag at not too high temperatures, T ^ 
uj c , is governed by the diffusive rectification! 81 1 which can be calculated quasiclassically using the local 
approximation for the density dependence of the conductivity. As a result, the sign of the drag at T ~ A 
oscillates but is opposite to what we found above for the ballistic regime: the drag is negative for equal 
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densities. (ST) At the "slopes" of the peak, p^ x scales with T and B in the following way 

( - T 2 ln(T£ 3 / 2 ), T«A, 

Pxx K { , , . (55) 

[ -T- Y B Z I 2 In 5, T> A, 

where the sign corresponds to the case of matching densities. 

Weakly ballistic regime, 1 <§; R c /a <C w c /A. This regime is qualitatively similar to the diffusive regime 
with 

r« A, 

A < T < TU = uj c {a/R c ) (56) 
T>co c (a/i? c ), 

The sign of the peak oscillates just like in the diffusive regime. 

Ballistic regime, lu c /A -C R c /a <C NA/u> c . The ballistic regime is most relevant experimentally. 
In this regime, the drag is governed by the particle-hole asymmetric effect of Landau quantization of the 
density of states and the sign of the drag oscillates, 



^xx 




Pxx <* 



T 2 B\n{B*/B), r«A, 

T -3 B 7/2 i n (£*/£), A « T « T, e A ln 1/2 (i? c A/aw c ), (57) 

_ T -i B 5/2 j T>T», 



where £?* ~ {mc/ e){v F / a 2 ^) 1 ^ (experimentally, the logarithmic factor in T* is typically of the order 
of unity, so that the intermediate regime may not be fully developed). We emphasize that the drag at low 
temperatures is positive for matched and negative for mismatched densities. 

Ultra-ballistic regime, NA/ui c <C R c /a. The drag for all temperatures is determined by the conven- 
tional contribution related to the curvature of the electron dispersion and is always positive, 



XX 




T< A, 

(58) 

T> A, 



At high temperature, T 3> u> c , the drag is governed by the conventional contribution (and is therefore 
positive) in all the regimes. It is linear in T in the diffusive regime (pf x oc TB- 1 ' 2 ). In all the ballistic 
regimes the drag resistivity scales as p® x oc T 2 B X / 2 for uj c <C T <C Vp/a and p® x oc TB 1 / 2 for 
T > v F /a. 

A comparison of Fig. [18] with Fig. 3 of Ref. ||79l reveals a remarkable agreement between the exper- 
imental findings and the theoretical results. In both the theory and the experiment, (i) p xx (T) shows a 
sharp peak at low temperatures; (ii) the sign of the drag in this temperature range oscillates as a function 
of the filling factor of one layer (at fixed filling factor of the other layer); (iii) the low-T drag is positive for 
equal filling factors and negative when the Fermi energy in one layer is in the upper half and in the other 
layer in the lower half of the Landau band; (iv) the high-T drag is always positive, independently of the 
difference in filling factors of two layers and increases monotonically with increasing T. Furthermore, it 
was observed by Muraki et al that in the low-temperature regime of initial increase of p xx , as well as in the 
high-temperature regime of "normal" drag, the drag resistivity can be described by an empirical scaling 
law, p xx oc (nj B)~ 2 7 f (T / B) . Theoretical results for both the low- and high-temperature regimes are in 
a nice correspondence with this prediction, with f(x) ~ x 2 . 
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Fig. 18 Schematic temperature dependence of drag in the ballistic regime 
for matched and mismatched densities. In the latter case the mismatch is 
chosen such that the drag is negative at low T (see text). Scaling of p® x 
with temperature in different regions is indicated: T <C A (I); A<T< 
T* (II); T, < T < u c (III), and T»iii c (IV). 

7 Photoconductivity 

7. 1 Microwave-induced magnetoresistance oscillations and zero-resistance states 

Recently, a number of new remarkable effects, important for both basic and applied physics, have been 
discovered in two-dimensional electron systems driven out of equilibrium by strong AC and DC fields. 
It was observed ll83l that dc resistivity p xx of a high-mobility 2DEG subjected to microwave radiation 
of frequency u> exhibits magnetooscillations with a period in ui set by the resonances with multiples of 
the cyclotron frequency uj c . Subsequent work on samples with an exceptionally high mobility has shown 
||84l[85l that for a sufficiently high radiation power the minima of these microwave-induced resistance os- 
cillations (MIRO) evolve into "zero resistance states"(ZRS), in which the dissipative resistance of a sample 
becomes vanishingly small. Unlike oscillatory p xx , the Hall resistivity p xy remained practically linear in 
u> c . A hallmark of these experimental findings is that the prominent oscillations of the photoconductivity 
dph — Pxx /p xy are observed at magnetic fields as low as 10 mT, and at relatively high temperatures up to 
~ 1 K, at which the Shubnikov-de Haas oscillations are completely suppressed. 

Presenting a novel class of magnetooscillations which lead, with increasing uj c , to apparently dissipa- 
tionless transport, the experimental results [83, 84, 85] have attracted much theoretical interest. In partic- 
ular, an explanation of the MIRO has been proposed 1 86 1 in terms of a combined effect of radiation and 
Landau quantization on elementary scattering acts for electrons colliding with impurities (in fact, a closely 
related theory was put forward long ago (87)). A systematic theoretical study of this mechanism of the 
MIRO (referred to as a "displacement" mechanism in what follows) was carried out in |88|. 

On the other hand, it was emphasized 1 89 1 that whenever the linear dc response theory predicts a nega- 
tive resistivity, this signifies an instability leading to the formation of domains of counter-flowing currents. 
The break-up of an ac-driven sample in current domains provides an explanation to the experimentally 
observed ZRS. 

A different mechanism of the MIRO, called here the "inelastic" mechanism, was proposed in lETTl and 
studied in more detail in 11901 |9T| (similar ideas were also discussed in |92|). The inelastic mechanism 
is associated with a radiation-induced non-equilibrium part of the distribution function of electrons f(e) 
which oscillates with varying e ± hw due to the Landau quantization. This mechanism yields the am- 
plitude of oscillations of the linear (with respect to the dc field) photoconductivity which is proportional 
to inelastic scattering time T m . The inelastic contribution dominates over the displacement one for r m 
larger than single-particle relaxation time r q , the condition which is fulfilled in the experiments. Apart 
from the magnitude of the effect, the two contributions are qualitatively different in their dependence on 
T and polarization of the radiation. In accord with the experiments, the inelastic contribution decreases as 
Ti n oc T~ 2 with increasing T and does not depend on the direction of linear polarization of the microwave 
field. By contrast, the displacement mechanism ll86l [87l [881 yields a T independent contribution which 
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depends essentially on the relative orientation of the microwave and dc fields, which clearly contradicts the 
experimental findings. 

7.2 Inelastic mechanism of MIRO 

We consider a high-mobility 2DEG, with r q <C rt r , subjected to a classically strong transverse magnetic 
field, u! c T tl - ^S> 1 (we use notations of Sec. [5] in particular, r q and r tr are specified in Eq. OTb . The 
photoconductivity <7 pn determines the longitudinal current flowing in response to a dc electric field £a c , 
j ■ £dc = fph^dc m tne P resence of a microwave electric field coscot. The more frequently measured 
ll83l l84l l85l [92 1 longitudinal resistivity, p pri , is given by p p h ~ pjLtfph, where p xy ~ eB/n e c is the Hall 
resistivity, affected only weakly by the radiation. 

Here we study the leading inelastic mechanism of MIRO thus taking into account only effects that are 
due to anon-trivial energy dependence of the non-equlibirum distribution function f(e). Photoconductivity 
Cph is given by the dc responce in the state with non-equlibirum f(e). According to Eq. (l38l l. 

Oph = -o£/<fei> 2 (e)&/(e), (59) 

where v(e) — v(e)/vq. The non-equlibirum distribution function f(e) is found as a solution of the sta- 
tionary kinetic equation 

£ l ^Ir5>( £±w ) [/(e±")-/(e)]+ £l 



'dc 9 



(60) 



On the right-hand side of Eq. (160) , inelastic processes are included in the relaxation time approximation 
(which is proven [91] to be sufficient under experimental conditions), and /t(e) is the Fermi distribution. 
The left-hand side is due to the electron collisions with impurities in the presence of the external electric 
fields. The first term describes the absorption and emission of microwave quanta; the rate of these transi- 
tions is proportional [61] tOT t ~g(e±o;) oc v(e±uj), see Eqs. (l33l and (f37T>- This term can be also extracted 
from the kinetic equation of Ref. [88 1. The second term describes the effect of the dc field and can be ob- 
tained from the first one by taking the limit to — > 0. Equation (|60T > suggests convenient dimensionless units 
for the strength of the ac and dc fields: 

r in /e£ u Df\ 2 ^c+^ 2 rs 2r in f eE Ac v F \ 2 f tt\ 2 

Ttr V U) J {U1 A —L0c) r tr \ U c J \UJcJ 

Note that and Qdc are proportional to r in and are infinite in the absence of inelastic relaxation processes. 

To first order in V u and Qdc — > 0, Eq. d60l l produces a non-equilibrium correction to /t(e)> /(e) — 
/t(e) = 0.25"P W 2^± ^ w ) ± w) — /t(e) ], which oscillates both with e/uj c and u>/w c due to 

£/w c -oscillations in the DOS. In turn, the oscillatory f(e) leads to tj/ti> c -oscillations of <j pri , Eq. (l59l . 

<W<^c = {v 2 {e))e+ (ujV u /A){P 2 (e)d E [D(e + cj)-P(e-uj)}} e (62) 

Here we took into account that the SdHO in the experiments are suppressed by temperature, T ^> To, so 
that, analogous to Eq. (|36"1 >. the energy integration results in averaging over e within the period oj c , denoted 
by the angular brackets. For separated LLs, cj c r q ^> 1, with the semielliptical DOS (l39l , Eq. d62l > gives 

ffph/og, = (16 Wc /3vr 2 r) {1 - VU^c/T 2 ) (w/r - nu c /T) + O {u c V u /T)]} , (63) 

47r$(a;) = a;Re[3arccos(|a;| - 1) - (1 + \x\)^/\x\(2- \x\)}. 

In the limit of overlapping LLs, the DOS is given by v = 1 — 25 cos with <5 = exp(— 7r/w c T q ) -C 1. 
The existence of a small parameter 5 allows one to calculate er p h to all orders in V u and Qdc, 

^^sin^+4Q d c 



^- = 1 + 2<5 2 

CT dc 



1 - 



1 + Vu, sin 2 2H + Q d 



(64) 
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Results (|64| > and (|63T > are shown in Figs. [T9landl20lfor several values of Vd, = V u \ Wc —q. 
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Fig. 19 Photoresistivity (normal- 
ized to the Drude value) for overlap- 
ping LLs vs Lo c juj at fixed coT q — 
2tt, and for different vL 0) = 
{0.24, 0.8, 2.4}. I - V characteris- 
tics at the marked minima are shown 
in Fig.HU 




Fig. 20 Photoresistivity (normal- 
ized to the Drude value) for sepa- 
rated Landau levels vs u) c /ui at fixed 
ajTq = 16-7T. The curves correspond 
to different levels of microwave 
power Pi 0) = {0.01, 0.03, 0.05}. 




Fig. 21 Current-voltage charac- 
teristics [ dimensionless current 
jx = (o-ph/crjjfdc vs dimen- 
sionless field £& c = Q dc ] at 
the points of minima marked by 
the circles in Fig.Qj?] The arrows 
show the dc field £ dc in sponta- 
neously formed domains. 



7.3 Zero-resistance states, strong-field domains, oscillatory compressibility 

At Vuj exceeding certain threshold value V*, er p h around minima becomes negative (see Figs. [19] and 
l20b . According to Ref. |89], this signifies an instability leading to the formation of domains with strong 
Hall field £d c = ±£dc anc ' counter-flowing currents. The value £| c is determined by the equation 
<J p h(Vu , £|c) = 0, which m the case of overlapping LLs, Eq. d64l ). has the solution (marked by arrows 
in Fig. ED 

'Pf n ] 1/2 ( a&™ ■ 2lTUJ ■ 2 7TUJ \~ 1 
— - 1 , V u = 4d — sin sin — . (65) 

V* \ V w c U c ) 

Equation d65l l relates the electric field formed in the domains (measurable by local voltage probe 1 93 ] ) with 
the excess power of microwave radiation. 

In the case of separated LLs, it suffices to keep the linear-in-T^ term only even for the microwave 
power > V* ~ T 2 /uui c at which the linear-response resistance becomes negative: The second order 
correction at V u ~ "P* is still small, UcV^/T ~ T/u> <C 1. The strength of the field in domains is 
determined by the scale at which inter-LL elastic scattering becomes efficient, £^ c ~ (u)\ / 'ev f) \J V tr / V q . 



Fig. 22 The microwave-induced correction to the compressibility 
(solid line) of a 2DEG as a function of u q /lj at fixed cuT q = 2-7T and 
microwave power V\ U( ,=o = 1. In the zero resistance state (ZRS), 
the electric field inside domains ££ c fixes the compressibility at the 
level shown by a dashed line. Inside the domain wall, the electric 
field £dc is smaller than ££ c and the compressibility depends on the 
local field as shown in the inset (for uj c /lj — 0.45, this ratio is 
indicated by the arrow. 



\/2 7r ev F 
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In Ref. ||94]| it was shown that the local compressibility of an irradiated 2DEG, \ = vq + S\, exhibits 
oscillations similar to the MIRO (see Fig.l22li. Calculation using Eqs. d60b and d34b yields 

f „V U — sin ^ +4Q dc 

S x /v = deD(e)d e [f T (e) - /(e)] = -S 2 "« , a "« . (66) 

J 1 + sin — + Qdc 

The key features of the effect are: (i) the period and the phase of the uj / oj c -oscillations in \ are the same 
as in o-ph, Eq. (1641) : (ii) the amplitude of the oscillations in x an d Cph have the same dependence on the 
electron temperature and microwave power; (iii) the ZRS corresponds to a plateau in the compressibility: 
inside the domains x = Vq{1— 25 2 )/2. Local measurements of the compressibility may provide a real space 
snapshot of the domain structure in the ZRS. Experimental work in this direction is currently underway. 

7.4 High-power effects, subleading mechanisms, fractional MIRO 
In addition to the peak-valley structure near integer lo/lo c (integer MIRO), several experiments l83l [85l 
|92l|95]|96l|97l|9U reported similar features near certain fractional values, uj/lu c = 1/2, 3/2, 5/2, 2/3.. 
("fractional MIRO", or FMIRO), which at elevated microwave power also evolved into ZRS ("fractional 
ZRS") |97l . Initially, FMIRO were ascribed to multiphoton processes [95 99]. This expanation, however, 
failed to reproduce the observations |98l , where the FMIRO only occured at ui below a certain threshold 
value. It was shown that the threshold can be explained in the framework of the single-photon inelastic 
mechanism 198111001 . Here, the FMIRO near combined resonances nui — mus c occured due to a resonant 
series of n single-photon transitions with real absorption (emission) of the microwave quanta [ 1 00 1 , in 
distinction to the virtual multiphoton processes. A systematic theory of the FMIRO II 1 II has shown 
that the existing theories 11991 1 1001 miss several important contributions. In particular, in the limit of 
well separated LLs the FMIRO are dominated by multiphoton inelastic mechanism. Provided Ti n /r q ^> 
1, the multiphoton displacement mechanism 1 99 1 yields a parametric ally smaller contribution and can 
be neglected. At weaker magnetic field the effects related to microwave-induced sidebands in the DOS 
become important. Close to the magnetic field at which the LLs start to overlap, the FMIRO are dominated 
by the single -photon inelastic mechanism [100|. Finally, in the regime of strongly overlapping LLs the 
FMIRO get exponentially suppressed. 

A unified picture of the photoresponce in the limit of overlapping LLs was recently developed in 111021 . 
On top of nonlinear interplay between the inelastic and displacement mechanisms at elevated microwave 
power, two novel mechanisms leading to the MIRO, "quadrupole" and "photovoltaic", were identified. In 
the quadrupole mechanism, the microwave radiation leads to excitation of the second angular harmonic 
of the distribution function. The dc response in the resulting nonequilibrium state yields an oscillatory 
contribution to the Hall part of the photoconductivity tensor which violates Onsager symmetry. In the 
photovoltaic mechanism, a combined action of the microwave and dc fields produces non-zero temporal 
harmonics of the stationary distribution function. The ac response in this state contributes to both the 
longitudinal and Hall MIRO. Provided Ti n /r q 3> 1, the inelastic mechanism still gives the dominant con- 
tribution to the diagonal part of the photoconductivity tensor. However, the quadrupole and photovoltaic 
mechanisms are the only ones yielding oscillatory corrections to the Hall part. Further, it was shown that 
a competition between various nonlinear effects (the feedback effects, the excitation of high angular and 
temporal harmonics of the distribution function, and the multiphoton effects) drives the system through 
four different nonlinear regimes with increasing microwave power. Most dramatic changes in the photore- 
sponse are due to the feedback effects. At V u ^> 1, the feedback from the microwave-induced oscillations 
of the isotropic part of the distribution /(e) leads to the saturation of the inelastic contribution, and to the 
strong interplay of the inelastic effect and all other contributions to the MIRO. In particular, the strong os- 
cillations of /(e) change sign of the most relevant parts of the displacement and photovoltaic contributions. 
At higher power, V u ^> Ld c T m , the feedback suppresses the effects on higher temporal and angular harmon- 
ics of the distribution function. At still higher power, V u Ti n /r q , the multiphoton excitation becomes 
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pronounced and starts to compete with the feedback effects. Finally, at V u 3> Tm/uj^T^, the feedback and 
multiphoton effects destroy all quantum contributions, restoring the classical Drude conductivity. 

7.5 New developments and open questions 

In spite of the essential advances in the understanding of nonequilibrium magnetotransport phenomena in 
a 2DEG, a number of questions remain open. A puzzling insensitivity of the MIRO to the direction of 
circular polarization of the microwave field was reported in [103]. The strong interplay [104] between 
the dc- 111051 H06l and microwave-induced oscillations deserves theoretical study. A further challenging 
direction of future experimental and theoretical research is the transition to the ZRS, the domain structure, 
electron transport and noise in the ZRS; first steps in this direction have been made in HI 071 110811 . In 
particular, the theory has not explained the seemingly activated temperature dependence of the residual 
resistance in the ZRS. Recently discovered £>-periodic magnetooscillations [ 1 09 1 which were ascribed to 
the microwave excitation of 2D edge magnetoplasmons warrant further investigation; in particular, their 
microscopic mechanism is still unclear. 
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